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' Abstract 

o . _ 

. We discuss the implications of the recent measurement of the Bs — Bs oscillation frequency AMs on the 

parameter space of R-parity violating supersymmetry. For completeness, we also discuss the bounds coming 
I from leptonic, semileptonic, and nonleptonic B decay modes, and point out some possibly interesting channels 

^ ! at LHC. 
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The Bs — Bs mass difference, recently measured by the DO [1] and the CDF [2] collaborations, is given by. 



O . 17 ps-i < AM, < 21 ps~^ {DO), AMs = {HMtoH ± 0.07) ps~\CDF). (1) 

, This result is consistent with the Standard Model (SM) prediction, which is estimated as 21.3 ± 2.6 ps^^ by 

r*^ ' the UTfit group [3] and as 20.9^^;^ ps'^ by the CKMfitter group [4]. The implications of AMs measurements 

Q-i, on the parameter space of New Physics (NP) have already been considered [5, 6, 7, 8]. However, given the 

^ • hadronic uncertainties in the SM prediction, along with additional uncertainties when NP is included, the 

O . present measurement of AMs does not provide a really strong constraint on NP [5, 6]. 

There have been some attempts to put bounds on the parameter space of R-parity conserving supersymmetry 
(RPC SUSY) from the AM, data [9]. In this paper, we would like to put bounds on the R-parity violating 
(RPV) SUSY couplings. We will use not only the AM, data but also the data on the leptonic, semileptonic, and 
nonleptonic branching ratios (BR) and CP asymmetries of B and Bs mesons that are affected by the particular 
RPV couplings. Such a work on B^ mesons may be found in [10] and this is an extension of that work to the 
Bs sector. For the relevant formulae, we refer the reader to [10]. 

It has been shown in [11] that RPV couplings involving sleptons (A and A' type) generate nonzero neutrino 
mass and one can put stringent constraints on them from the WMAP data [12]. The exact bounds depend on 
the relation of the mass matrices with the CKM matrix. Anyway, such a study forces us to consider only those 
couplings which can still be relatively large, and in this paper we derive better bounds on some of these product 
couplings than those coming from [11]. 

A major motivation for this study is the Bs physics that is going to be probed at LHC-b, and even at CMS 
or ATLAS during the low-luminosity run of the Large Hadron Collider (LHC). The leptonic and semileptonic 
decays are clean and any enhancement over the SM expectations will signify some NP. Also, the phase x 
Bs — Bs mixing comes in the subleading order of the CKM matrix and is expected to be very small (~ 0.03), 
so any CP-violating effect in Bs decays not involving an up quark will be interesting. 

Effects of RPV SUSY on B physics have been discussed extensively in the literature [13, 14, 15, 16, 17]. 
Constraints coming from neutral meson mixing have been discussed in [10, 18, 19]. However, in these papers. 
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the Bg sector could not be dealt with, since only the lower bound on AM^ existed then. This paper, in a sense, 
is the completion of the series. All the computational details that have been taken into account in [10, 19] {e.g., 
the NLO QCD corrections for short-distance effects, inclusion of both SM and RPV) are incorporated in this 
paper. 

The paper is arranged as follows. In Section 2 we outline the relevant formulae necessary for the analysis, and 
give the numerical inputs in Section 3. The analysis on Bg — Bg box and the decay processes is in Section 4, 
and we conclude and summarize in Section 5. 



2 Basic inputs 
2.1 Bg — Bg mixing 

The off-diagonal element M12 in the 2x2 effective Hamiltonian causes the Bg — Bg mixing. The mass difference 
between the two mass eigenstates AM^ is given by (following the convention of [20]) 

AM, = 2|Mi2|, (2) 

with the approximation \Mi2\ 3> |ri2| (this seems a good approximation even for the Bg system). If we have n 
number of NP amplitudes with weak phases one can write the mass difference between mass eigenstates as 

AM, = 2[|Mf2^|2 + ^ |Ml2|2 + 2|Mf2^| ^ \Ml^\ cos2{0sm -0i) + 2j2J2 I^i2ll^l2l cos2(^,- - 9^)]'^^ (3) 

i i i j>i 



For Bs — Bg system, the short-distance SM amphtude is 

MSm ^ = ^^{y,gV:,fris^rnBjlBB.ml,So{x,). (4) 

where generically Xj = rrtj/m^, fs^ is the Bg meson decay constant, and 773^ and Bb^ parametrize the short- 
and the long-distance QCD corrections, respectively. The function So is given by 

^ , , Ax-llx'^+x^ 3x^lna: , , 

= 4(1 -x)^ - 20^- 

If the NP amplitude has a nonzero phase, then there will be an effective phase in Bg — Bg mixing amplitude, 
whose presence may be tested in the hadronic B factories. In the presence of NP, the general AF = 2 effective 
Hamiltonian can be written as 

5 3 

'^efr' = E ci(M)Oi(M) + E ci(M)Oi(M) + H.c. (6) 

i=l i=l 

where /i is the regularization scale. The effective operators Oi and Oi are given in [10]. The Wilson coefficients 
c, at = include NP effects, coming from couplings and internal propagators. However, for most of the 
NP models, and certainly for the case we are discussing here, all NP particles are heavier than mw and hence 
the running of the coefficients between mw and ^ = 0{m},) are controlled by the SM Hamiltonian alone. In 
other words, NP determines only the boundary conditions of the renormalization group (RG) equations. For 
the evolution of these coefficients down to the low-energy scale, we follow Ref. [21], which uses, for B^ — B'^ 
mixing, /z = nib = 4.6 GeV. The expectation values of these operators between Bg and Bg at the scale fj, are 
analogous to those as given in [10]. The Bb^{h) parameters have been taken from [22]. 
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2.2 R-parity violation 



It is well known that in order to avoid rapid proton decay one cannot have both lepton number and baryon 

number violating RPV couplings, and we shall work with a lepton number violating model. This leads to both 
slepton (charged and neutral) and squark mediated decays, and new amplitudes for — Bg mixing, where 
particles flowing inside the box can be (i) charged sleptons and up- type quarks, (ii) sneutrino and down type 
quarks, (iii) squarks and leptons. One or both of the scalar particles inside the box can be replaced by W 
bosons, charged Higgs bosons and Goldstone bosons (in a non-unitary gauge) (see Fig. 1). We follow the usual 
practice of avoiding the so-called "pure SUSY" contributions to the box amplitudes, i.e., those coming from 
charginos, neutralinos or gluinos inside the loop. Not only the strongly interacting superparticles are expected 
to be heavier than the electrowcak ones (and hence the contribution being suppressed), but also one can choose 
SUSY models where these contributions become negligible {e.g., alignment in the squark sector, or Higgsino- 
dominated lighter chargino, to kill off the respective boxes.) Since the current lower bound on the slepton mass 
is generally weaker than that on squark mass by a factor 2-3, the slepton mediated boxes have greater chance 
to be numerically significant. 

We start with the superpotential 

Wa' = KokLiQiDl, (7) 

where i, j, A: = 1, 2, 3 are quark and lepton generation indices; L and Q are the S'C/(2)-doublet lepton and quark 
superfields and is the SU (2)-singlet down-type quark superfield respectively. Written in terms of component 
fields, this superpotential generates six terms, plus their hermitian conjugates: 

Adndi + 44^1 + (4)*^14 - eldnui - uid^L - (4)*e>il + H.c. (8) 



-LQD — ^ijk 



With such a term, one can have two diff'erent kind of boxes, shown in Fig. 1, that contribute to Bg — Bg mixing: 
first, the one where one has two sfermions flowing inside the loop, alongwith two SM fermions [23], and secondly, 
the one where one slepton, one W (or charged Higgs or Goldstone) and two up-type quarks complete the loop 
[18]. It is obvious that the first amplitude is proportional to the product of four A' type couplings, and the 
second to the product of two A' type couplings times Gf- We call them L4 and L2 boxes, respectively, for 
brevity, where L is a shorthand for A'. 

We will constrain only products of two A'- type couplings at a time, and assume a hierarchical structure, i.e., only 

one product is, for all practical piirpose, simultaneously nonzero (but can have a nontrivial phase). This may not 
be physically the most appealing scenario but keeps the discussion free from unnecessary complications. This 
product can in general be complex. Any product is bounded by the product of the bounds on the individual 
terms, which we call the direct product bound (DPB). Interesting bounds are those which are numerically 
smaller, and hence stronger, than the corresponding DPBs. The DPBs are mostly taken from [11], and we 
highlight those products which are more tightly constrained than their respective DPBs. The detailed formulae 
of the box amplitudes may be found in [10]. 



2.3 Semileptonic and leptonic decay channels 

The RPV couplings that may contribute to Eg — Eg mixing should also affect various E decay modes. Let us 
first consider the leptonic and semileptonic modes. 

The expected BRs of leptonic flavour conserving A_B = 1 processes within SM are much below the experimental 
numbers (except B K^*H'^i~), so one can safely ignore the SM effects as well as the R-conserving SUSY 
effects to put bounds on the RPV couplings. (The flnal state leptons must be the same if the product coupling 
contributes to Bg — Bg mixing.) The leptonic decay modes are theoretically clean and free from any hadronic 
uncertainties. The semileptonic modes have the usual form-factor uncertainties. 

To construct four-fermion operators from A' type couplings that mediate such semileptonic and leptonic B 
decays, one needs to integrate out the squark or slepton field. The product RPV coupling may in general 
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Figure 1: R-parity violating contributions to Bg — Bs mixing. Figure (a) corresponds to L4, while figure (b) to 
L2 amplitudes (see text for their meanings). For L4-, there are similar diagrams with squarks and leptons (both 
charged and neutral), as well as diagrams with left-chiral quarks as external legs and quarks and sneutrinos 
flowing in the box. For L2, there are diagrams where the W is replaced by the charged Higgs or the charged 
Goldstone. The internal slepton can be of any generation, and so can be the internal charge +2/3 quarks, 
generically depicted as u. 

be complex. However, since all the leptonic decays are one-amplitude processes (only RPV, for all practical 
purposes) there is no scope for CP-violation; one can only look at nonzero BRs. By the same argument, we can 
take all couplings to be real without any loss of generality. 

The effective Hamiltonian is of the form [24] 

-\C3klm [Pj-f^PLl'l] \d^li^PLdm] + B.C., (9) 

where 

3 \/* \/ 3 \i* y 

We take any one to be nonzero at a time, mi, is the left-chiral up/down squark mass, taken to be degenerate. 

The entire leptonic B^ decay amplitude is solely due to new physics, as far as detectability is concerned. In 
RPV models, squark-mediated A' A' type interactions are responsible for such purely leptonic decays. As already 
pointed out, the bounds are robust in a sense that they are free from any theoretical uncertainties (except for 
the decay constants of i?s), and do not depend on the phase of the RPV couplings. For the B mesons, no such 
leptonic mode has yet been observed. The corresponding upper limits on the BRs are of the order of 10~^ 
iox I = \x and 10~^ ior i = e modes. With 300 GeV squarks, from the bounds that one obtains here, a BR 
at most of the order of 10^^ can be expected. Thus, we do not expect to see such leptonic channels before 
the next-generation hadronic machine or super e+e~ B factories. However, a number of semileptonic modes 
b — > s£~^i~ have been observed and the BRs are at the SM ballpark. 

The decay width of Bg — > £~£'^ is given in [24]. For the semileptonic decays, we use 
{K{p2)\srb\B{p,)) 
(.^(P2,e)| V^tA^ \BM) 



= P^F,{q^) + q^ "'^ (Fo(g^) - Fi(g^)) , 

= \ [-tl^(9')£/..a/3e*"'P"/ 
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±Aoiq^)iP-q)e;±A+{q'){c*-p,)P^ 

±A.{q^){e*-pM (11) 

where ms, and are the meson masses, Pi{p2) is the momentum of the initial (final) meson, e is the polar- 
ization vector of the vector meson 0, P = Pi + P2, Q = Pi — P2, = q^ltiqi, A^ = q^lp-l^li^ ^o,± and 
P(i/o){(f) are the form factors. The values of these form factor are taken from [25, 26]. 

The RPV matrix elements for the decay mode -B(s/d) — * {4>/K)l'^l~ (where I = e,jjb): 

M(i?(,/d) ^ {ct>IK)lttr) = iB,3,2 il+^i^PLr) (((</./i^)|(;7d)7^Pfl6|i?(,/rf))) (12) 



2.4 Nonleptonic decay channels: 

There are four types of slepton-mediated nonleptonic decays that can proceed through the relevant RPV cou- 
plings. Among them, h ccs and b sss rates are bound to be undetectably small if these couplings are to 
be compatible with the neutrino mass bounds [11]. We, therefore, focus on the b uus and b dds type 
transitions. They mediate the channels B — > ttK, B pK. The corresponding Bg decay channels do not have 
data at comparable level. Note that both BRs and CP asymmetries for these channels have been measured 
[27, 28]. While the data does not uniquely point to NP, the trend is encouraging. 

We will use the Conventional Factorization (CF) model [26] to analyze the effect of RPV SUSY on the ttK and 
pK channels. While the validity of such a simple approach may be questioned, it is not wildly off the truth, at 
least for these channels. The effective Hamiltonian reads 

Hrpv = dfk^{d„i^PLdj)sidk-i^,PRb)s + dfi^^idr,j''PLb)s(dkl^PRdj)s 

+uf^k{^nl''PLUj)i(dkl^PRb)s (13) 



where 



jR _ ^ijk^ in3 ,L _ ^»3fc^ inj R _ m3 



Following the standard practice we shall assume that the RPV couplings are hierarchical i.e., only one combina- 
tion of the couplings is numerically significant. Lot us assume, to start with, only d\i2 and ufi2 to be nonzero. 
The QCD corrections are easy to implement: the short-distance QCD corrections enhance the {S— P)y. {S + P) 
RPV operator by approximately a factor of 2 while running from the slepton mass scale (assumed to be at 100 
GeV) to rub [29]. 

The RPV amplitude for B — > ttK is given by 

M^oK- = -L[«f,,(-i?,(AW)+42^i-) + 4,,i-(42^)] (15) 
= ^[4i2(-i?i(4'l) + 4'!^) + ^12^(4'^)] (16) 

M^.j^ = d^,,{R^A^^^) (17) 

=<i2(-^i4'^) (18) 
The expressions for the B — !■ pK amplitudes will be similar to those shown above, with Ri replaced by i?2 and 
replaced by We use the shorthand 

Rl=2-. ; 77 r, -^2 = 2-^ ; r, (19) 

[mu + ms){mb - mu) (m„ -h ms)(m6 - to„) 
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and 



a'-;^ = 2fKmX^'{rnl){ep.pK), A% = 2/,m,Ff ^^(m2)(e,.p^). (20) 



3 Numerical Inputs 



Quantity 


Value 


AM, 


Z 1 i\ ' j ' J r ^ 

(17.31+^;^^ ± 0.07)ps-^ 


7 


50-72 




0.55 ±0.01 




166 GeV 




4.23 GeV 


nibimb) 


4.6 GeV 


mc{mb) 


1.3 GeV 


nidirrib) 


5.4 McV 


ms(2 GeV) 


125 MeV 


fB\/BB,\jLQCD 


(0.245 ± 0.02i;° °°3) GeV 


IK.I X IQi 


9 979+O.dl 
''-0.01 


\Vos\ X IQi 


9.73 


l^tsl X 103 


41 61+°-i2 


IKbl X 103 


4.4 ±0.3 


|T46| X 103 


42.0 ±0.7 


IVibI X lOi 


9.99 



Table 1: Input parameters used for the numerical analysis, from [2, 4, 21, 30]. 

The major sources of the numerical inputs are: (i) the Heavy Flavor Averaging Group (HFAG) website [31] for 
the latest (summer 2006) updates on B physics; (ii) Particle Data Group 2006 edition [32]; and (iii) the inputs 
used in the CKMfitter package [4]. The quark masses and Wilson coefficients have been taken from [21, 33]. 
We use the following numbers. 

The masses for all the mesons B^, B~, w, p, and K are the corresponding central values as given in [32]. The 
meson decay constants (in GeV) are: 

U = 0.133, fK = 0.158, fp = 0.210. (21) 
The transition formfactors [34] at = are given by 

Fo{B ^ K) = 0.38; Fo{B -> tt) = 0.33; Ao{B ^ p) = 0.28, (22) 

and ^b(O) = i^i(O). 

The quark masses have been evaluated in the MS scheme. The pole mass for the top quark is about 5 GeV 
higher and the mass for the bottom quark is 4.6 GeV. The CKM elements arc shown in Table 1. 

The leptonic and semileptonic BRs for the B meson, which are of interest to us, are as follows [1, 32, 31]: 

Br(B Kl+r) < (0.57 ± 0.07) x 10"^ {I = e/p); 
Br(S^Xe+e~) = (0.55 ± 0.09) x 10"^; 
Bv{B^Kp+p-) = (0.61 ± 0.08) X 10"^; 
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Bt{Bs e+e") 



< 4.1 X 10"®; 

< 1 X 10""^; 

< 5.4 X 10"^; 



For the ttK and pKmodes, the data reads [32, 31, 35, 36, 37]: 

Br{B° ■K-K+) = (18.9 ± 0.7) x 10"^ 

Br(B° ^ 7r°if°) = (11.5 ± 1.0) X 10"*^ 

Br{B+^7T+K") = (24.1 ± 1.3) X lO"'^ 

Br{B+^Tr°K+) = (12.1 ± 0.8) x IQ-*^ 

Br{B° ^ p-K+) = {9.9t\i) x IQ-*^ 

Br{B° p°K°) = (5.1 ± 1.6) x IQ-^ 

Br{B+ ^ p°K+) = (4.231^:^?) x 10"'^ 

A'^%{B° ^n^K"^) = 0.115 ±0.018. (23) 



To evaluate the QCD corrections, we take as{m\) = 0.1187 [32], and take the SUSY scale Ms = 500 GeV. 
The precise value of this scale is not important, however, and we can take it to be at the squark mass scale 
(300 GcV) without affecting the final results. The exact cvohition matrix can be found in [21] and [33]; for 
our purpose, it is sufficient to note that for the Bg system, the operator Oi is multiplicatively renormalized by 
a factor 0.820 at the scale p, = 2 GeV, and the operator O4 at mw changes to (2.83O4 + O.O8O5). We again 
stress that theoretically the procedure is questionable for boxes with light quarks flowing in the loop. However, 
the numbers that we obtain are fairly robust and one can very well drop the NLO corrections altogether, if 
necessary, without compromising the results. Since the O5 admixture is small, one can take the central values 
for these parameters without introducing too much error. The relevant S-parameters for Bg system are taken 
from [22] 

We take all sleptons to be degenerate at 100 GeV, and all squarks at 300 GeV. We also take tan/3(= V2/V1) =5 
(very low values are excluded by LEP, and the numbers are not sensitive to the precise choice of tan/3), and 
the charged Higgs boson mass as 200 GeV (lower values are disfavored from b — > 57). 



4 Analysis 

4.1 Bg — Bg mixing 

For the B^ system, the bounds are summarized in Table 2. When the product coupling is complex, we show only 
the real part, since the bound on the imaginary part is almost equal to this. The reason is easy to understand: 
the bounds are obtained when the RPV coupling has a phase opposite to that of the SM coupling, so that the 
interference is destructive. At the limit where the RPV coupling determines the mixing amplitude, the phase 
is irrelevant. The effect of this destructive interference is clear in Fig. 2(a) and Fig. 2(b). For a more detailed 
explanation, we refer the reader to [10]. 

The relative magnitudes of the bounds are also easy to miderstand. For example let us consider the bounds 
on A A vis-a-vis A -22 A -33. The relevant box diagrams have the same particle content; but the first one is 
proportional to VtsVcb (~ 0{\'^)), and the second one to VtbVcs (~1)- The relative suppression in A enhances 
the limit on the RPV coupling. 

Though most of the bounds are of the same order in magnitude, these are, theoretically, an improvement over 

those obtained earlier [18, 23, 19]. We have taken into account all possible amplitudes (and the interference pat- 
terns play a nontrivial role), including the SM one, but have systematically neglected the pure supersymmetric 
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A'V 


Only 


Complex, 


combination 


real 


real part 


(132) (i33) 


1.01 X K)-^ 


1.0 X 10"^ 


(i22)(i23) 


(S.2 X K.)--' 


8.0 X 10"'^ 


(il2)(il3) 


1.2 X 10"^ 


3.7 X 10-^ 


(i22)(i33) 


1.8 X 10"^ 


1.7 X 10-1 


(132) (123) 


3.2 X 10"* 


3.0 X 10"* 


(i22)(il3) 


3.47 X 10-"^ 


3.45 X 10"^ 


(il2)(i23) 


5.16 X 10"^ 


7.56 X 10-^ 


(i32)(il3) 


1.4 X 10"^ 


1.3 X 10"^ 


(il2)(i33) 


9.0 X 10-1 


9.0 X 10-1 


(123) (133) 


1.66 X 10"^ 


5.1 X 10"^ 


(122) (132) 


1.66 X 10"^ 


5.1 X 10"^ 


(i21)(i31) 


l.()() X 


•1.1 X 



Table 2: Bounds on A' A' combinations from Bg — Bg mixing. The table displays the magnitudes only, and not 
the signs. 




-0.01 0.01 

Re[(i12)(i13)l 




(a) 



(b) 



Figure 2: (a) Allowed parameter space for \'i\2^'ii3 (t>) The allowed paramater space of RPV phase for 
|A^i2'^ji3l'''^^ich can gives CP assymetries in the Bg — > K-k decay. 
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A'A' 
combination 


Quark 
level 


Meson 
level 


Bound 


(i21)(i31) 


b — > dds 


Bd TT^RO 
B- 'K°K- 


2.45 X 10-3 


(il2)(il3) 


b dds 

b — > uus 


B- ■k'^K- 

Bd n+K~ 


2.82 X 10-3 


(i22)(i32) 


b sss 


Bs -> <pK, 


2.33 X 10-3 



Table 3: Some of the possible nonleptonic transitions mediated by the RPV couplings discussed in the paper. 



boxes coming from gaugino exchange. The reason is that those boxes decouple in the heavy squark limit, and 
one can always take an RPV model embedded in a minimal supersymmctric theory where such FCNC processes 
are somehow forbidden. It was shown in [19] that the bounds are fairly robust even if one takes into account 
such SUSY contributions. Furthermore, the QCD corrections are implemented upto NLO. 

4.2 Nonleptonic decay channels 

As we have mentioned, the B ttK and B pK numbers are encouraging for NP enthusiasts. However, the 
hadronic uncertainties are significant. Also, one must have a nonzero strong phase between the SM tree and 
the SM penguin amplitudes to explain the direct CP asymmetry data on B — > Tr+if-. Thus, it is of importance 
to explore the data in conjunction with Bg — Bg mixing. In our analysis, we use the CF model, as mentioned, 
and vary the strong phase difference between the SM tree and the RPV amplitudes from to 2tt. To take 
into account the hadronic uncertainties, wc (i) vary the SM amplitude from its CF value by 20%, and (ii) vary 
the RPV weak phase in the range [0:27r]. The allowed parameter space of RPV couplings, in the magnitude- 
phase plane, comes out with separate island-like structures. For example, A^3]^A^2i li^s between 2.05 x 10^^ to 
2.45 X 10^3 and the corresponding phase lies between 104° to 120°, whereas A-j3A-;^2 li^s between 2.31 x 10-^ 
to 2.82 X 10-3 and the corresponding phase lies between 85° to 105°. We quoted the highest value in Table 3. 
These bounds are much stronger then that coming from Bs — Bg mixing. 

We have not considered the B (77, 'r]')K modes. Though they are mediated by 6 ^ sqq (q = u, d) transitions, 
the BRs for those modes cannot be explained simultaneously by such a simple new physics structure [38]. 
Similarly, the mode B (f>K* has not been considered, since the longitudinal polarization anomaly for this 
mode cannot be explained without a contribution from tensor current but RPV SUSY do not provide for such 
tensor current structures, at least at the tree-level [39]. However, we note that unless a product coupling is at 
least of the order of 10-3, the RPV contribution is unlikely to affect the SM amplitude. 

4.3 Leptonic and semileptonic decay channels 

The detail procedure of the leptonic and semileptonic decay channels are given in [24] . While their bounds were 
for 100 GeV squarks, we scale the numbers for 300 GeV squarks. The bounds are shown in Table 4. Note that 
the channel B Kfi'^iJ.~ gives the best bound. Also, there is no bound involving rs in the final state, but 
we can estimate the number of Bg t^t~ decays. The relevant RPV coupling is AJj^jA'^^^. It can easily be 
checked that unless i = 1, the product is so constrained from neutrino mass [11] that even at LHC-b, there is 
no hope to detect a RPV signal in this channel. For i = 1, the coupling A312A313 should be less than 2.8 x 10-3, 
which in turn translates into a bound on the BR to be less than 2.7 x IQ-^. Note that the SM expectation is 
about 7 X 10-'^. 
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A'A' 
combination 


Quark 
level 


Meson 
level 


Bound 


(2i2)(2i3) 


b — > s/i+/i~ 


Bs jj+pT 


7.7 X 10"^ 


(2i2)(2i3) 


b sfi'^^~ 


Bs (/i/i+Zi" 


4.9 X 10-^ 


(2i2)(2i3) 


b S/Lt+/Z~ 


B KjjL+n- 


6.6 X 10-"* 


(Ii2)(li3) 


b — > se+e" 


B Ke+e- 


7.7 X 10"* 



Table 4: Some of the possible leptonic and semileptonic transitions mediated by the RPV coupling relevent 
with Bs — Bs mixing are given here. It has shown that in many cases the bounds coming from this decays are 
much better then coming from mixing. 



A'A' 
combination 


Related 
process 


Current 
bound 


Previous 
bound 


(112) (113) 


Bs Ke+e- 


7.74 X 10"* 


1.52 X 10-^^ 


(122)(123) 


Bs Ke+e- 


7.74 X 10-'* 


9.7 X 10-«^ 


(132)(133) 


Bs -> Ke+e- 


7.74 X 10-"* 


6.0 X 10-^ 


(212)(213) 


Bs Kfi+fi' 


6.57 X 10"* 


1.52 X 10-^ 


(222)(223) 


Bs K^i+n- 


7.74 X 10-4 


9.7 X 10-«^ 


(232)(233) 


Bs ^ Kfi+fi- 


7.74 X 10-'* 


6.0 X 10-",, 


(312)(313) 


B- TT^K- 
Bd ^ n+K- 


2.8 X 10-'^ 


1.52 X 10-*^ 


(322)(323) 


Bs — Bs 


8.2 X 10-3 


9.7 X 10-6^ 


(332)(333) 


Bs — Bs 


1.01 X 10-2 


6.0 X 10-^ 


(122) (133) 


Bs — Bs 


1.8 X IQ-i 


3.75 X 10-9^ 


(i32)(i23)t 


Bs — Bs 


3.2 X 10-4 


1.56 X 10-*^ 


(i22)(il3) 


Bs — Bs 


3.47 X 10-2 


9.75 X lO-'^^ 


(il2)(i23) 


Bs — Bs 


5.16 X 10-2 


1.52 X 10-^^ 


(i32)(il3) 


Bs — Bs 


1.4 X 10-3 


1.56 X 10-*^ 


(il2)(i33) 


Bs — Bs 


9.0 X 10-* 


5.8 X 10-^ 


(123) (133) 


Bs - Bs 


1.G6 X 10-2 


5.8 X 10-^, 


(i22)(i32) 


Bs oKs 


2.3 X U)--' 


1.0 X 10-'^ 


(121) (131) 


Bd - n^K" 
B- Tr^K~ 


2.45 X 10-3 


1.56 X 10-^^ 



Table 5: Bounds on real A'A' combinations from Bs — Bs mixing and correlated decay channels. The DPBs, 

displayed in the last column, occur from neutrino constraints with no mixing scenario [11]. The product marked 
with a dagger is bounded from tree-level Bs — Bs mixing(~ 0(10-'')). 
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4.4 Comparison between bounds coming from mixing and decay 

In Table 5 we summarize our results, displaying the bounds on all A' A' type products that may be responsible 

for Bs — Bs mixing. We find that a number of them may have better bounds from semileptonic or nonleptonic 
decay modes. While the neutrino constraints are indeed tight, we obtain a tighter constraint for most of the 
products. 



5 Summary and Conclusions 

In this paper wc have computed the bounds on the product couplings of the type A'A' coming from Bg — Bg 
mixing. Though such an analysis is not new, we have implemented several features in the analysis which have 
not been taken into account in earlier studies. Previously there was a lower limit on AM;,, here wc have used 
the current bound on it and considered the exact expression for the box amplitude taking all possible processes, 
including that from SM. The QCD corrections to the amplitudes have been taken upto the NLO level. We 
have considered the possibility that the RPV product couplings may be complex. The analysis is done in the 
benchmark point m.H+ = 200 GcV, tan/3 = 5, all sleptons degenerate at 100 GeV and all squarks degenerate 
at 300 GcV, and n(\Q,icc;ting the pure; MSSM contribution to the box amplitudes (by possibly applying some 
underlying FCNC suppression principle, like alignment of the squark mass matrices). 

It is to be observed that in some cases, our bounds are actually weaker than those obtained earlier by saturating 
the mass difference with RPV alone. The reason is that destructive interference with the SM amplitude plays 

a very crucial role in determining the bounds, particularly when the phase of the RPV coupling is arbitrary. 
There is an intricate interplay among different amplitudes as can be seen in Fig. 2. 

In some cases the bounds obtained from semileptonic B decays are better. Of course, one can enhance the 
squark mass to a limit where these bounds become weaker then those obtained from the box (the latter is not 
much affected by decoupling the squarks), but such extremely massive squarks are not interesting, even for the 
LHC. 

However, some of these couplings may affect the nonleptonic decay modes (which, being slcpton mediated, 
cannot be suppressed by decoupling the squarks). For some cases, the bounds coming from such decays are 
tighter than those coming from mixing. 
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